Wave packet dynamics with Bose-Einstein condensates 
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We study wave packet dynamics of a Bose condensate in a 
periodically shaken trap. Dichotomy, that is, dynamic split- 
ting of the condensate, and dynamic stabilization are analyzed 
in analogy with similar phenomena in the domain of atoms in 
strong laser fields. 
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Recently it has become possible to prepare Bosc- 
Einstcin condensates of alkali eases [| with up to 10 
magnetically trapped atoms. The condensed state is a 
macroscopically populated quantum state well localized 
in the magnetic trap. It is, therefore, an ideal tool to 
study wave packet dynamics under experimentally fea- 
sible conditions. There are many interesting quantum 
phenomena resulting from the electronic wave packet dy- 
namics; in this Letter we refer, in particular, to the phe- 
nomena of wave packet dichotomy and stabilization ex- 
hibited by an electron bound by an atomic potential in 
presence of a strong laser field We argue that the 
same phenomena occur in the dynamics of the conden- 
sate wave function. The analogy is based on the fact that 
in the frame of reference moving with the free electron 
oscillating in the field, the Kramers-Henneberger frame, 
the effect of the laser is equivalent to a periodic shak- 
ing of the atomic potential along the laser polarization 
axis. A condensate in a periodically shaken trap could, 
therefore, a priori show a similar behavior. 

As the intense laser field drives the electron, the pro- 
cess of ionization of the atom occurs. By increasing the 
laser intensity, one normally increases the ionization rate. 
However, for very intense fields of high frequency, this 
rate eventually starts to decrease with intensity - this is 
called atomic stabilization 0]. In this process the elec- 
tronic wave packet remains bounded, i.e. well localized 
in space (without spreading), although highly distorted 
due to the combined effects of the laser field and the 
atomic potential . This effective atom-laser potential ex- 
hibits a double well structure which splits the electronic 
wave packet into two spatially separate parts; this effect 
is called dichotomy. To achieve stabilization it is neces- 
sary to turn on the laser adiabatically in order to ensure 
that the atomic ground state will evolve to the ground 



state of this effective atom-laser potential. This type of 
stabilization |3j has never been observed experimentally 
P|, since it requires very intense high frequency fields, 
which currently can only be generated in a form of a 
very short pulse; as an electron in an atom is highly un- 
stable, it would thus be most likely ionized during the 
turn-on of such a pulse. 

Let us analyze the analogy between the electron and 
the condensate in more detail. The electron bound by an 
atomic potential U(r) and interacting with a laser field 
of amplitude £e z (polarized along the z-direction) is, for 
our purposes, best described in the Kramers-Henneberger 
frame of reference, in which the interaction with the laser 
field results in an effective time dependent "atomic" po- 
tential: 
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+ U e (r + a L sm(u) L t)e z ) 



V B (?,t) = 0, (1) 



where oll = (e£) / '(m e w^) is the electron excursion am- 
plitude, while ul is the laser frequency. 

Consider now a condensate with N atoms in a mag- 
netic trap V{f) which is periodically shaken along the 
z-axis. In a Hartree-Fock treatment, the state of the 
condensate is described by the Gross-Pitaevskii equation 
(GPE) ||. It accurately describes the wave function of 
the condensate "J in presence of particle interactions in 
thermal equilibrium at temperatures well below the criti- 
cal temperature. Furthermore, the time dependent GPE 
describes the dynamics of the condensate in more general 
time dependent conditions J^-Q]. We, therefore, have: 
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+ V(r + a(t)e z )+gN\V(r,t)\ 2 



*(f,i) = 0. 
(2) 



Here a(t) = ao sin(wt) is the shaking amplitude. If we 
identify the amplitude ao and the frequency u> of the 
shaking with the electron excursion amplitude and 
the laser frequency ojl respectively, Eq. (|2|) is very similar 
to Eq. ([!]). An important difference is the presence of a 
nonlinear coupling term with a coupling constant g — 
Airh 2 a s /m, where a s > is the s-wave scattering length. 
The effects due to the presence of the nonlinear term 
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as well as the larger mass of the atom will be discussed 
below. We model the trapping potential by 



(3) 



for V(r) < V c , and V{r) = V c otherwise. This is a har- 
monic potential with frequencies Q x ,y,z which is cut at an 
energy V c . It thus resembles U e (r), in particular we have 
that V(f — > oo) like U e (f — > oo) is constant. This trap- 
ping potential can be realized as an effective (adiabatic) 
trap potential "dressed" by the microwave coupling be- 
tween a trapped and an untrapped state. We will discuss 
this model later on in the paper. Another possible real- 
ization would be a condensate in a dipole trap potential 
formed by a strong off-resonant laser held. The condensa- 
tion might take place in such a trap, or the magnetically 
trapped condensate may be loaded into it. Note that in 
absence of the cut V c , i.e. for an exactly harmonic poten- 
tial the shaking of the trap would lead to an undistorted 
oscillation of the condensate wave function. 

If the time scale l/u>z, in Eq.(jl]) is shorter than the 
typical time scale of the electronic motion, the combined 
laser-atom potential can be replaced by its time aver- 
age over one laser period; ideally the wave function of 
the electron will then evolve adiabatically from a state in 
the bare atomic potential into the corresponding state of 
this time-averaged atom-laser potential For a large 
enough amplitude oil, this time average potential will 
have a double well structure and thus the electron wave 
function will exhibit dichotomy. Using the same reason- 
ing for the condensate we replace the GPE (Eq.(||)) by 
a GPE with a time-averaged potential if the time scale 
~ 1/ui of the shaking is shorter than all other relevant 
time scales: 
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V eS (r,a )+gN\y(r,t)\ 2 



#(f,t) = 0, 
(4) 



where the time-averaged potential is given by 
1 



V c g(f,a ) 



2n 



dipV(r + «o sva(ip)e z ); ip = Lot. (5) 



For sufficiently large shaking amplitudes, a double well 
structure appears as illustrated in Fig. 1(a). We therefore 
expect dichotomy as well as stabilization to appear also, 
due to the trapping of the atomic wave function in the 
two potential wells of V c s near the turning points of the 
oscillation. Indeed, for large enough ao, stable states are 
dynamically formed in the double well. Even though the 
particles in such states are sometimes out of the trap, no 
"ionization" occurs. This is due to the fact that the time 
scale for particle motion is much larger than the time 
scale l/u> of shaking, so that the particles do not have 
enough time to react on being momentarily out of the 
trap. Due to the large atomic mass, the requirements on 
uj are in fact much less stringent than in the electronic 
case. 



The initial state of the condensate for temperature 
T ~ is given by the lowest eigenstate \l/o of the time 
independent GPE 



fl/V 2 
2m 



+ V(r) + Ng\Hf (f)\ 2 ]V (f) = 0. (6) 



where /i is the chemical potential. By replacing V by V e e 
in Eq. (||) we determine the ground state of the system in 
presence of the double well potential. The corresponding 
chemical potential is /i e g. 
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FIG. 1. (a) Cut harmonic potential(solid line) and time av- 
eraged potential (dashed line) for ao = 30a z and V c = 80?ifi z ; 
a z = (h/2mS~l z ) 1 / 2 and gN = 100. (b) Dressed state poten- 
tials (V+ and V-) and time averaged potential of the lower 
one (V-) from Eq.(|); a = 30a z , lu r = 100fi z , A = 200fi 2 . 

In contrast to Eq. ([!]), the GPE accounts for atomic 
interactions, so that the condensate splitting depends on 
gN. With increasing N the wave function will at some 
point overcome the potential barrier between the wells in 
V e fi, and the splitting will disappear. This will happen 
when the effective chemical potential (/i c ff ) in the double 
well potential V e g exceeds the height of the double well. 
Therefore the number N is limited from above. 

The validity of the above picture depends on the turn- 
on time of the shaking. Ideally, an adiabatic switching 
will transfer the ground state of the GPE (Eq.(||)) at 
t = to the ground state of the GPE with V replaced 
V c ff. Too rapid passage from the cut-harmonic potential 
into the double well form induces transitions to the higher 
states of V e g and even into the continuum. To avoid these 
transitions, we turn on the shaking gradually, setting 
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a sin 2 ( f tM sin(wi), for < t < t on 
a(t)=i - (7) 



ao sin(a;i), 



for t > t n 



The turn-on time t on must be 3> 2tt/uj; in practice 
we take t on > (50 — 150) x 2tt/lo ~ 27rao/a z ti;, where 
a 2 = \JXij1m£l z . 

A typical experimental sequence to demonstrate the 
above phenomena will consist in the following. Initially, 
the system is prepared in the state v£/ of Eq.(||). At 
time t = we slowly start shaking the trap along the z- 
axis (Eq.(^)). Let us first assume that the perpendicular 
motion does not play a significant role, and consider the 
motion along the z-axis only. In Fig. 2(a) we show the 
the time evolution of the ID-condensate wave function. 
The parameters used are the same as for Fig. 1(a), i.e. 
lo = 10ft z , V c = 80hn z , a = 30a z , and gN = 100. 
The gradual splitting of the condensate wave function 
can be clearly seen. Moreover almost all particles remain 
trapped. 

Due to the cut V c in the trapping potential (Eq.(||)) 
atoms may escape from the trap. However, as can bee 
seen from Fig. 2(a), almost no atoms do it; this indi- 
cates stabilization. To study this phenomenon in more 
detail we lower the cut-off energy V c to 50hil z to favor 
the escape from the trap. The escape rate is calculated 
using standard absorbing boundary conditions (b) ; we 
observe a decrease in the escape rate if the shaking ampli- 
tude ao, increases. More specifically, increasing ao from 
15a 2 to 20a z , the escape rate decreases by a factor ~ 2. 
Due to the large atomic mass the escape rate is over- 
all very small, i.e. < 1% of the trapped population per 
100 shaking cycles. For the same reason, the condensate 
stabilization occurs already for relatively small shaking 
frequencies huj < V c — /i c ff ■ The condensate stabilization 
is, therefore, more pronounced and can be achieved in ex- 
perimentally more accessible conditions than its analog 
in strong field dynamics of electrons. 





FIG. 2. (a) Time evolution of the condensate density 
^(z,*)! 2 undergoing 300 shaking cycles; ao = 30a z , 
uj = 10f2 z , t on = 150 x (27r/w), V c = 80hn z . The 
nonlinear coupling is gN — 100, that corresponds to 
hfj, — 14.13Q Z ; (b) Same as (a) for the total conden- 
sate density |*(z,t)| 2 = \^!(z,t)\ 2 + \^ (z,t)\ 2 , calculated 



from the ID two state model; ao = 30a z 
ton = 150 x {2-k/u), ujr = 100f2 z , A = 200ft z . 
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We turn now to a discussion of a model with two inter- 
nal states, a trapped state (F, mp ^ 0) and an untrapped 
state (F'^mp = 0) . They are coupled via a microwave 
field which allows coherent transitions between the states 
[^j . In the rotating wave approximation the GPE is given 
by 
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. 5 A(|*o(r,i)| 2 + |*i(r,i)| 2 ) 



mn 2 z (z + a{t)) 2 /2 hu R /2 
Hujr/2 hA 



*o(M 
*i(r,t) 



0. (8) 



Here ^0,^1 are the wave functions of atoms in the 
trapped and untrapped state normalized to the respec- 
tive fraction of atoms in these states, A is the detuning of 
the microwave from the transition frequency, and ljr is 
the Rabi frequency of the microwave transition. For sim- 
plicity we assume all coupling constants equal to g. For a 
large enough Rabi frequency ljr, the coupled states can 
be replaced by uncoupled dressed states with energies: 



V±( z ,t) = - [mn 2 z (z + a(t)) 2 /2 + hA 
± \J (mn 2 z {z + a(t)) 2 /2 - hAf + h 2 u 



2 , ,2 
R 



(9) 



The potential V_ resembles our model potential; this 
is illustrated in Fig. 1(b) in which we also plot time 
averages of the time dependent dressed state potential. 
Obviously Eq.(||), with V replaced by V_, is not exact, 
since it neglects entirely the non-adiabatic transitions to 
an upper branch of the dressed potential. The model is 
nevertheless reasonable, since the exact numerical treat- 
ment based of the two-component GPE leads to very 
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similar results, as shown in Fig. 2(b). In the model of 
Eq. (||) we account rigorously for non-adiabatic (Landau- 
Zener) transitions from the lower to the upper "dressed" 
state manifold. Obviously, some of the atoms may cease 
to oscillate in the lower potential and instead they may 
be transferred to the upper "dressed" potential, destroy- 
ing the dichotomy. In order to avoid the Landau-Zcncr 
transitions w has to be smaller than a certain critical 
value H . On the other hand for the validity of the time- 
averaged potential uj has a lower limit of approximately 
2f2 z . In consequence dichotomy and stabilization in a 
realistic system occur in a limited range of w's. 

Finally, we have generalized our study to a 3D case. 
We chose parameters that resemble those of the MIT 
experiment [0 , that is a cigar shaped trap with a small £l z 
and equal frequencies perpendicular to it. We assume the 
shaking occurs along the long z-axis. As shown in Fig. 3 
the presence of perpendicular motion does not invalidate 
the conclusions from the ID approximation: the splitting 
of the condensate is clearly visible. 
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FIG. 3. 2D cross-section of the condensate density 
\^(x,z,t)\ 2 for a cigar shaped trap at i = 400(27r/u;); 
a = 60a z , u = 10n z , t on = 250 x (27r/w), V c = 30htt z ; 
gN = 100, and Q x = Q y = 5Q Z . 

One should mention that the completeness of our anal- 
ysis requires to check the stability of the solutions of the 
time dependent GPE. Unstable solutions might lead to 
a depletion of the condensate |l(| . In a first attempt we 
verified that the solutions of an oscillating condensate in 
a harmonic trap lead to stable solutions of the GPE, and 
therefore to no depletion. 

Summarizing, we have shown that Bose-Einstein con- 
densates are ideal tools to study wave packet behavior 
first predicted in the realm of atoms in superstrong laser 
fields. We have shown that both, dichotomy and stabi- 
lization against atom escape can be achieved in conden- 
states. Although these effects are very difficult to realize 
in electron-atom systems, their observation in conden- 
sates seems quite feasible. Furthermore, we believe that 
wave packet dynamics of condensates might lead to inter- 



esting possibilities of condensate state engineering. For 
instance, so far double peaked condensates have been cre- 
ated using laser "knives" that cut a single condensate into 
two parts. We offer here an alternative method to achieve 
a similar dichotomy in a more controlled way which opens 
new perspectives for condensate interference studies. 
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